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EVALUATION OF ANALYTICAL PROCEDURES FOR PREDICTION 
OF TURBULENT BOUNDARY LAYERS ON A POROUS WALL 


by Charles E. Towne 
Lewis Research Center 

SUMMARY 

An analytical study has been made to determine how well current boundary layer 
prediction techniques work when there is mass transfer normal to the wall. The data 
that were considered in this investigation were for two-dimensional, incompressible, 
turbulent boundary layers with suction and blowing. Some of the bleed data were taken 
in an adverse pressure gradient. An integral prediction method was used with three dif- 
ferent porous wall skin friction relations, in addition to a solid-surface relation, for the 
suction cases. A numerical prediction method was also used. Comparisons were made 
between theoretical and experimental skin friction coefficients, displacement and mo- 
mentum thicknesses, and velocity profiles. The integral method with one of the porous 
wall skin friction laws gave very good agreement with data for most of the cases con- 
sidered. The use of the solid -surface skin friction law caused the integral method to 
overpredict the effectiveness of the bleed. The numerical technique also worked well 
for most of the cases. 


INTRODUCTION 

There are many situations where turbulent boundary layer prediction techniques 
must account for mass transfer normal to the wall. One of these is transpiration cool- 
ing, which involves blowing into the boundary layer. Another is boundary layer control 
in adverse pressure gradients by means of suction. This latter case is particularly ap- 
plicable to flow in supersonic inlets, where bleed is often necessary to prevent separa- 
tion (refs. 1 and 2). It is not difficult to include the terms involving mass flux normal to 
the wall in the boundary layer equations, but little attempt has been made to systemat- 
ically compare how well current prediction techniques work under these conditions. 



This report, therefore, presents the results of two analytical procedures, along 
with data, for several turbulent boundary layers on a porous wall. The integral method 
is that of Sasman and Cresci (ref. 3), modified to include the terms involving the veloc- 
ity normal to the wall. With integral procedures, a modification to the usual solid- 
surface skin friction relation to account directly for nonzero velocities at the wall is also 
necessary. Three different skin friction laws (refs. 4 to 6), applicable to boundary lay- 
ers with suction, are used in the integral analysis. The numerical method of McDonald 
and Fish (ref. 7) is used to verify the integral results and is also used with some injec- 
tion cases. Comparisons are made between theoretical and experimental displacement 
thickness, momentum thickness, and skin friction coefficient. Velocity profiles are also 
presented for the numerical method. 

Analytical methods require knowledge of the mass flow distribution normal to the 
wall through the bleed or blowing region. Because of a lack of data of this type for com- 
pressible flow, the data considered in this report are for incompressible, two- 
dimensional flow. Several different mass transfer rates and distributions are investi- 
gated. Some of the data that are used were taken in a pressure gradient. 


SYMBOLS 


A function of Re^ (see eq. (2)) 

B function of Re^ (see eq. (2)) 

Cj skin friction coefficient 

C, skin friction coefficient with v = 0 
to w 

H shape factor, 6*/ 6 

h enthalpy 

j dimensionality parameter 

"wake” value of mixing length 
M Mach number 

Pr laminar Prandtl number 

p pressure 

R radius of axisymmetric body 

Re x Reynolds number based on x 

Re Reynolds number based on 0 

U 

T temperature 
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u 

u 

V 
v 
X 
x 

Y 

y 

z 

r 

A 

6 

6 * 

6 


transformed x -velocity 
x -velocity 

transformed y -velocity 
y -velocity 

transformed x -coordinate 
coordinate along body surface 
transformed y-coordinate 
coordinate normal to body surface 

UlRe e 

ratio of specific heats 
transformed boundary layer thickness 
boundary layer thickness 


displacement thickness 


•/ 


1 - -£UL- ]dy 
^ Pe u e/ 


momentum thickness 


, f 

L P e u e \ 


dy 


pi viscosity coefficient 

v kinematic viscosity coefficient 
p density 

r shear stress 

Subscripts: 

e local external conditions 

i transformed value 

s local stagnation conditions 

w conditions at wall 

0 free stream stagnation conditions 

Superscript: 


evaluated at reference conditions (see eq. (A5)) 



METHOD OF ANALYSIS 


Two analytical techniques, one numerical and one integral, were used in this inves- 
tigation. The integral method that was used is essentially that of Sasman and Cresci 
(ref. 3). This method uses the transformed momentum and moment of momentum inte- 
gral equations, requiring starting values of momentum and displacement thicknesses. 
The equations are coupled and solved simultaneously by a Runge-Kutta integration rou- 
tine. For this study the equations were modified to include terms involving the velocity 
normal to the wall and to allow the use of different skin friction relations. The details 
of this modification are given in the appendix. 

The skin friction relation used in reference 3 is 


C.= 0.246 e- 1 - 561H i( Vi 

1 V V n 


■0. 268 T . x 0.268 

T0 


(1) 


Using this solid-surface relation in a prediction method with suction will lead to inaccu- 
rate predictions of the boundary layer growth, as pointed out in reference 4. It is there- 
fore necessary to use a skin friction equation that directly accounts for the velocity nor- 
mal to the wall. 

Three different relations of this type were tried. The first was from Thompson 
(ref. 4). It is based on a three -parameter family of velocity profiles, and is applicable 
to incompressible flows with small adverse pressure gradients. This relation is pre- 
sented both graphically and as an analytic approximation by Krishnamurthy in ref- 
erence 4. The analytic approximation was used in this study, and is of the form 


C f = e AH+B (2) 

Letting z = In Re fl , for v /u = 0 the values of A and B are given by 

(7 W C 

A = 0.019521 - z[0. 386768 - z(0. 028345 - 0. 0007017z)], 

B = 0.191511 - z[0. 834891 - z(0. 062588 - 0.001953z)] 


For v w /u g = -0. 005 

A = 1. 07085 - z[0. 831747 - z(0. 106843 - 0. 004428z)] , 
B = -2.338049 + z[0. 015834 - z(0. 047968 - 0.003308z)] 
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For v w /u e = -0.010 


A = -1. 00747 - e (1 °- 97531 - 8 - 1080555 In z) 

B = 110.92 e " 106z -2.94 

The value of Cj at a given v w /u g was found by linear interpolation or extrapolation 
using these values. 

The other two skin friction relations that were used are based on incompressible 
flat plate data. Both depend on Cj 0 , the value of the skin friction coefficient without 
bleed. One of these equations (ref. 5) is 



-1.62 


/3 w v w 

Pe u e 



1,62 p w v w \ 
c fo Pe u e/ 


The other, from reference 6, is 


C 


f 



-0.748 


e 


p w v w 1 
Pe u e C fo 


For these equations the value of Cj Q was found by using the Sasman-Cresci method with 
equation (1). 

The numerical method is that of McDonald and Fish (ref. 7). It is based on a mixing 
length turbulence model and uses a general implicit finite -difference scheme to solve the 
boundary layer equations. In reference 7 it is stated that the value of the outer layer 
normalized mixing length, l^J 6, is allowed to vary in the streamwise direction for non- 
equilibrium boundary layers, requiring the use of the integral turbulence kinetic energy 
equation. For this investigation, however, 1^/5 was set at a constant value of 0 . 09 . A 
correction to the turbulence model for low Reynolds number flows, given in reference 8, 
was included in the analysis. 

The numerical method requires starting velocity and temperature profiles. Several 
different techniques can be used to obtain these profiles. The experimental profiles can 
be used, if available, or they can be found by using the flat plate turbulent boundary 
layer profile family of Maise and McDonald (ref. 9). In addition, an "equilibrium” start 
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can be used in which an iteration upon the starting profiles is performed. This method 
tends to relax and improve a poor guess for the starting profiles, and was used for all 
cases in this study. 


RESULTS AND DISCUSSION 
Zero Pressure Gradient 

The zero pressure gradient data that were used are those of Simpson (ref. 10). 
These data were taken on a flat plate with a smooth porous surface of sintered bronze. 
The porous region was 45. 7 centimeters wide and 244 centimeters long. A boundary 
layer trip was located 0.318 centimeter upstream of the start of the porous surface. 
Boundary layer profiles were mea sirred along the plate centerline at various distances 
from the leading edge. Table I presents the free stream conditions and bleed rates for 
the Simpson cases used in this study. Some of the mass transfer rates are functions of 
x, the distance along the plate from the start of the porous region. The run numbers 
correspond to those of reference 10. 

The zero pressure gradient bleed cases were analyzed using the integral method of 
reference 3 with four different skin friction laws. The starting values of displacement 
and momentum thickness were those that gave good agreement at the first experimental 
probe location for the zero mass transfer case (run 31067). The theoretical values of 
the skin friction coefficient for run 31067 were used for Cj Q in the correlations of ref- 
erences 5 and 6. 

The results are shown in figure 1. The use of the solid-surface skin friction rela- 
tion of reference 3 tended to overpredict the bleed effectiveness. The skin friction laws 
of references 5 and 6, based on flat plate data, gave somewhat better results, but the 
best results were obtained by using the curve fit to the skin friction law of reference 4. 
An exception was the high bleed case shown with an expanded scale in figure 1(b). Here 
the experimental momentum and displacement thicknesses were close to their laminar 
asymptotic values of 0.012 and 0.025 centimeters (ref. 11). This may indicate that the 
boundary layer flow is becoming laminar due to the large bleed rate. Under these con- 
ditions it is difficult to say which skin friction relation gave the best agreement. When 
the correlation of reference 5 was used for this case, the integral thicknesses went to 
zero because the bleed terms in the governing equations became dominant. 

All 11 of the zero pressure gradient cases were analyzed using the numerical pre- 
diction technique. The starting value of displacement thickness for all cases was the 
same as that used with the integral method. The initial profiles were found from ref- 
erence 9 as previously discussed. 
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Figure 2 compares theoretical and experimental values of skin friction coefficient, 
momentum thickness, and displacement thickness. The configurations with uniform 
blowing and bleed are presented in figures 2(a) and (b), those with variable mass transfer 
rates in figures 2(c) and (d), and those with step changes in the mass transfer rate in 2(e) 
and (f). The results for the no mass transfer case are shown on each plot as a reference. 
Excellent agreement was obtained for all configurations. Velocity profiles are presented 
in figure 3 for five of these cases. Again excellent agreement was obtained. 


Nonzero Pressure Gradient 

The adverse pressure gradient bleed data were taken on the flat part of the upper 
surface of a symmetrical airfoil by Thompson (ref. 12). The airfoil had a chord length 
of 81.2 centimeters and a circular leading edge of radius 8.25 centimeters. The suction 
surface itself was 54.6 centimeters long and made of a smooth porous plastic material. 
Boundary layer profiles were measured at seven stations along the airfoil, with the first 
just upstream of the start of the suction surface. In addition, static pressure taps were 
located along the surface to determine the external velocity distribution. 

Five different boundary layers were analyzed, identified as layers B, D, F, G, and 
H in reference 12. Layers B and D were measured with the airfoil at an angle of attack 
of -6° for layer B and +5° for layer D. For layers F, G, and H, a leading edge section 
was attached, and the airfoil was used as a flap, at a flap angle of approximately 17°. 
Figure 4 shows the experimentally determined external velocity distributions and figure 5 
the bleed flow distributions for the five cases that were studied. These experimental 
values were those used for the analysis. 

The results of the integral analysis of reference 3 were compared with the adverse 
pressure gradient data. The experimental momentum and displacement thicknesses just 
upstream of the bleed region were used as starting values. The method was first used 
without bleed to determine the values of Cj 0 for the skin friction correlations of ref- 
erences 5 and 6. The results are shown in figure 6. Note that layers F, G, and H ap- 
pear close to separation without bleed. 

Figure 7 compares the results of the integral method using the various skin friction 
relations with the adverse pressure gradient data. For all cases, the solid -surface skin 
friction relation overpredicts the effectiveness of the bleed, and to a greater degree than 
for the zero pressure gradient cases previously discussed. The relation of reference 4 
generally gave very good agreement with the data. The correlations of references 5 
and 6 gave poor results. This should be expected, since these correlations are based on 
zero pressure gradient data. 
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For layer B, the momentum and displacement thicknesses, shown on an expanded 
scale in figure 7(a), went to zero because of the high bleed rate when the skin friction laws 
of references 3 to 5 were used. The correlation of reference 6 gave too high a value for 
C|, causing the method to blow up just downstream of the start of suction, and is not 
shown in figure 7(a). The correlations of references 5 and 6 both caused the method to 
eventually blow up for the remaining cases also. 

For the numerical technique, the experimental velocity profile just upstream of the 
bleed region was used to start the calculations. The starting temperature profile was 
determined by assuming a Crocco temperature distribution. The equilibrium starting 
method was used, as previously discussed. 

Figure 8 compares the results of the numerical technique with the experimental 
values of momentum and displacement thickness. The theoretical skin friction coeffi- 
cient is also shown, but no experimental skin friction measurements were made. Good 
agreement was obtained for the three lower speed cases, layers B, D, and H. The 
agreement was not as good, however, for layers F and G. The same results were found 
for the velocity profiles, figure 9. There are several possible reasons for the disagree- 
ment found in layers F and G. Both of these cases involve a fairly strong adverse pres- 
sure gradient and a small bleed rate. There may be limitations in the mixing length 
model that was used in this study when applied under these conditions. Small errors in 
the experimental bleed rate may be important in an adverse pressure gradient. There 
was also a blockage effect due to the supporting structure beneath the porous surface that 
was not taken into account. 

For a further evaluation of the numerical method of reference 7, two additional 
cases were considered, both involving injection in a pressure gradient. The data were 
taken by McQuaid (ref. 13) on a flat plate with a porous region 76. 2 centimeters long. 

One of the cases, called pressure distribution I in reference 13, was for a slightly ad- 
verse pressure gradient. The other, pressure distribution II, was a favorable gradient. 
Boundary layer measurements were made with a pitot pressure probe at various stations 
along the plate. Static pressure taps were used to determine the streamwise velocity 
distributions. 

These velocity distributions, along with the experimental injection rates, were used 
as input for the analysis. The values of the velocity were all of the order of 15 meters 
per second. The local injection rate v w / u e was around 0.002 for pressure distribu- 
tion I and 0.008 for pressure distribution II. The calculations were started at the first 
probe location by using the experimental displacement thickness. The boundary layer 
profiles were found by the method of reference 9, as previously discussed. 

Figure 10 presents the theoretical and experimental values of momentum and dis- 
placement thickness for both pressure distributions. The predicted skin friction coeffi- 
cient is also shown. The velocity profiles are given in figure 11. Excellent agreement 
was obtained for both cases. 
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SUMMARY OF RESULTS 


Comparisons of analysis and data have been made for turbulent boundary layers with 
suction and blowing normal to the wall. Two different prediction methods were used. In 
addition, three different porous wall skin friction relations were used with an integral 
method for comparisons with the suction cases. Results of a numerical technique are 
also presented. The data were for incompressible, two-dimensional flow with several 
different mass transfer rates and distributions. Pressure gradient effects were present 
in some of the data. Comparisons were made of theoretical and experimental skin fric- 
tion coefficients, displacement and momentum thicknesses, and velocity profiles. 

The results of the investigation can be summarized as follows: 

1. In the zero pressure gradient cases, the integral theory gave better agreement 
with data when the porous wall skin friction relations were used in place of the solid - 
surface relation. One of the porous wall equations, intended for use in zero and mild 
pressure gradients, was consistently better than the others. 

2. For the cases with pressure gradient, the two porous wall skin friction laws 
based on flat plate flows gave poor agreement with data. The third, applicable to cases 
with mild pressure gradient, gave good agreement for the low speed cases with high 
bleed rates , but slightly underestimated the displacement thickness for the cases with 
lower bleed rates. 

3. The numerical technique gave excellent agreement with data for all the zero pres- 
sure gradient cases considered. Good agreement was also obtained for the four low 
speed adverse pressure gradient cases. The velocity profile predictions, however, were 
not as good for the two more severe pressure gradient cases. 

In order to apply theory to boundary layers across bleed regions in supersonic in- 
lets, more experimental work needs to be done in compressible flow and in adverse pres- 
sure gradients. Skin friction data is needed in addition to velocity profiles. Detailed 
measurements of the bleed mass flow distribution are also required. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 6, 1974, 

501-24. 
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APPENDIX - MODIFICATION OF SASMAN-CRESCI INTEGRAL EQUATIONS 

For steady, compressible, turbulent flow the time averaged boundary layer equa- 
tions are 


3M + £(£v) + jPudR = o (Al) 

dx dy R dx 


p u** + pv 

3x 3y dx 3y 


(A2) 


where j = 0 for two-dimensional flow and 1 for axisymmetric flow. 

These equations are transformed using the method of reference 3, given by 


and 


where 



2 T, 


1 - =11^ + 0.22 Pr 1 / 3 + (0.5 - 0. 22 Pr 1 / 3 ) _® 


The continuity equation (Al) becomes 


au + av 
ax 3Y 


. U dR 
+ ] - — = 
R dX 


0 


and the momentum equation (A2) transforms to 

u3U + v SU = ^u dU e t 1 T 3t 
ax 3Y h 0 e dx p e T e 3Y 


(A3) 


(A4) 


(A5) 


(A6) 


(A7) 
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Using equations (A6) and (A7) and following the procedure outlined in references 14 
and 15, but without setting V w = 0, gives the transformed momentum integral equation 


d^ + V_U_ e 

dX U dX 
e 




l]dy| + j!L^-I^ = ^i-^L 

R dX U T T n „ tt 2 
e e 0 p e U e 


(A 8) 


and the transformed moment -of -momentum integral equation 


dH^ 

_ _1_ 

d^e 

+ ‘) 2 ®i - » 

dX 

U 

e 

dX 

2 




2(H. - 1) 
H?(H. + l)ef 



1 YdY 



Now following reference 3, these become 


d^ = 

e i 


,-i^SS + 

1 R dx 

dx 

" M e 

*■ K \ j 

and 





Pe /^e\ V2 % + /T s \ (y+1)/2(y ' 1) C, 
L p 0 T w\ T 0/ u e \ T 0/ 2 


(A 10) 
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dx 


J_^e 

2M g dx 


[ H t( H i + l) 2 (Hi - «] 


1 + 1 — - 11 


H f + 4H. - 1 


T q j (H. + 1)(H. + 3)J 

2 rr, /m \(r+l)/2(r-l) 


- Mil (Hj + 1 )(h? - l)M if— ^ 
*i Vl H i ) T\r o y 


K - >) 

hi 


r+l)/2(y-l) h 2 _ l 


V2 v 

W 

•i ' 

\ T oj 

2 + «i Po V 

K t o) 

' u e 


These equations are in a slightly different form than those of reference 3. The momen- 
tum integral equation has been left in terms of 0. , and l^/hg has been replaced by 
T w /T 0 . The skin friction coefficient appears explicitly in equations (A7) and (A8), al- 
lowing the use of different equations for Cj when v w is nonzero. The relation between 
v w and V w is given by 



T T 0 V w 


T T u 
e w e 


(A 12) 
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TABLE I. - FREE STREAM CONDITIONS AND BLEED RATES 


FOR SIMPSON DATA (REF. 10) 


Run number 

Velocity 

Temperature 

Pressure 

Bleed or injection rate, 


V 

m/sec 

T 

e’ 

K 

Pe ’2 

N/m z 

v /u 
w 7 e 

31 067 

EQ 

292. 2 

100 900 

0 

122 866 

US 

304. 1 

102 500 

+0.00099 

122 066 

14. 10 

304.7 

101 200 

+0.0038 

21 067 

12.96 

291.4 

101 600 

-0.0012 

2 867 

13.02 

291.0 

102 200 

-0.0046 

12 467 

13. 12 

289.9 

100 200 

-0.056 Re'^/u 

X G 

3 867 

12. 80 

291.6 

100 900 

+0.101 Re' 1 / 5 /*! 

x e 

61 367 

12. 83 

293.7 

100 500 

+5.0X10" 5 x 

42 567 

13. 18 

291. 8 

102 200 

r 0. 0 for Re x < 8. 14xl0 5 
\+0. 001 for Re x > 8. 14xl0 5 

5 167 

13.26 

292. 6 

102 000 

f 0.0 for Re x < 8. 15X10 5 
1+0.0039 for Re >8. 15X10 5 

62 367 

13. 12 

294. 6 

100 700 

f -0. 0024 for Re x <7.95xl0 5 
lo.Ofor Re x >7.95xl0 5 













Displacement thickness, Momentum thick- Skin friction 
6*, cm ness, 0, cm coefficient, C> 



(c) Run 12467; v w /u e = -0. 056 Re^ 5 /u e . (d) Run 62367; step in suction rate at x = 92 centi- 

meters from v w /u e = -a 0024 to 0. a 

Figure 1. - Concluded. 
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<r . 012 


Theory Data 
(ref. 7 ) (ref. 20) 


Run Vyy/Ug 








(a) For uniform blowing. 

Figure 2. - Theoretical and 





(b) For uniform suction. 


results for Simpson flat plate. 





Displacement thickness, 6, cm Momentum thickness, 0, cm Skin friction 

coefficient, C 







displacement thickness, 5*, cm Momentum thickness, 0, cm .Skin friction coefficient, 



20 
















Distance along surface, cm 
(e) Layer H. 

Figure 7. - Concluded. 





Figure 8. - Theoretical and experimental results for Thompson airfoil. 





Figure 8. - Continued. 




Distance along surface, cm 


\9 







m ‘aoejjns o* ibiujou aoueisio 


29 




Pressure 



Distance along surface, cm 

Figure 10. - Theoretical and experimental results for 
McQuaid pressure distributions I and II. 



Page Intentionally Left Blank 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON. D.C. 20546 


OFFICIAL BUSINESS 
PENALTY FOR PRIVATE USE *300 


SPECIAL FOURTH-CLASS RATE 
BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



POSTMASTER : 


If Undeliverable (Section 15S 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information co?icerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 





